Relativistic analysis of magnetoelectric crystals: 
extracting a new 4-dimensional P odd and T odd pseudoscalar from Cr 2 03 data 
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Earlier, the linear magnetoelectric effect of chromium sesquioxide Cr203 has been determined 
experimentally as a function of temperature. One measures the electric field-induced magnetization 
on Cr203 crystals or the magnetic field-induced polarization. From the magnetoelectric moduli of 
Cr2 03 we extract a 4-dimensional relativistic invariant pseudoscalar 5. It is temperature dependent 
and of the order of ~ 10 -4 Yq, with Yq as vacuum admittance. We show that the new pseudoscalar 
5 is odd under parity transformation and odd under time inversion. Moreover, 5 is for Cr203 what 
Tellegen's gyrator is for two port theory, the axion field for axion electrodynamics, and the PEMC 
(perfect electromagnetic conductor) for electrical engineering. 
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I. INTRODUCTION 

Our paper addresses the magnetoelectric (ME) effect. 
This effect has been established since the 1960's in Cr203 
crystals (see the reviews by O'Dell [1| and, more recently, 
by Fiebig @). The ME effect, in linear approximation, 
is described by magnetoelectric susceptibilities or moduli 
that have been measured by different groups. 

It has been a long-standing discussion whether these 
moduli fulfill a certain condition, as predicted by Post in 
1962 [3j . This condition was dubbed Post constraint by 
Lakhtakia, see, e.g., [3,[5j]. Numerous arguments against 
the validity of the Post constraint were put forward, some 
of them are mentioned in [|| and 0], e.g.. However, in the 
end, we must turn to the experiments and their proper 
evaluation. 

Following Post 0], we provide here a relativistic in- 
variant formalism of the electrodynamics of moving me- 
dia. The violation of the Post constraint is measured 
by a pseudoscalar modulus, which has not been deter- 
mined so far. If this pseudoscalar vanishes, the Post con- 
straint is fulfilled, otherwise it is violated. On the basis 
of experimental data, we determine this pseudocsalar (or 
axion) piece of the magenetoelectric moduli and find it 
non-vanishing. Therefore the Post constraint cannot be 



uphold as a general valid relation. 

However, our paper has also an interdisciplinary pur- 
pose. Our result of the non-vanishing pseudoscalar pro- 
vides a physical structure that also shows up in the the- 
ory of electric networks, more exactly in the theory of two 
ports, as Tellegen's gyrator Q, in electrical engineering 
as perfect electromagnetic conductor (PEMC) [9(, and 
in elementary particle physics as the hypothetical ax- 
ion field [ll[. These interrelationships support each 
other. Since the axion in elementary physics is the only 
left hypothetical object in this context, our results make 
also the existence of the axion particle more likely. 

In Sec. 2, we give a short description of the ME effect. 
In Sec. 3, a 4-dimensional electrodynamic framework for 
moving media is built up and the electromagnetic con- 
stitutive tensor introduced for local and linear media. In 
Sec. 4 we discuss Dzyaloshinskii's theory [l2j of 0^03 
and extract therefrom the mentioned pseudoscalar. In 
Sec. 5 we finally determine the pseudoscalar for the first 
time and discuss some of its properties in Sec. 6. Then, in 
Sec. 7, we turn to the interdisciplinary part and discuss 
the existence of the pseudoscalar for network theory, elec- 
trical engineering, and elementary particle physics. In 
the concluding section, we collect our results. 



II. MAGNETOELECTRIC EFFECT 
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In classical electrodynamics for a local linear medium, 
which is at rest in the reference frame considered, the con- 
stitutive law reads D = eeqEi and H = B/(/x/Xo)- Here £o 
is the electric constant (permittivity of free space) and 
Ho the magnetic constant (permeability of free space), 
whereas e and /i are the (relative) permittivity and per- 
meability, respectively, of the medium under consider- 
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ation. Further more, the admittance of free space is 
Y = = y^So/no, with flo as vacuum impedance, 

and the speed of light c = 1/y/eoVo- 

If an external B field in some suitable medium induces 
an electric excitation D and an external E field a mag- 
netic excitation H, the constitutive law mentioned has 
to be extended by so-called magnetoelectric pieces, see 
O'Dell The general local and linear constitutive law, 
if the medium is anisotropic, reads 



D = (e)e E+ (ai)Y Q B, 



(1) 
(2) 



We have to read (TT]) and ^ as tensor equations, with 
(e), (cki), and {0-2) as dimensionless 3x3 matri- 

ces. Hence we expect 36 permittivity, permeability, and 
magnetoelectric moduli in general. The constants £o> 
and fio are required for dimensional consistency. 

The existence of nonvanishing (cti ) and (012) matrices 
was foreseen by Landau-Lifshitz [13j for certain magnetic 
crystals and proposed by Dzvaloshinski i fl2l| specifically 
for the antiferromagnet Cr203. Astrov [lj] (for an elec- 
tric field) and Rado & Folen [15j (for a magnetic field) 
confirmed this theory experimentally for Cr2C>3 crystals. 
For reviews, see [l6j — there other magnetoelectric crys- 
tals are listed, too — and . 



III. FOUR-DIMENSIONAL MAXWELLIAN 
FRAMEWORK 



In order to extend the formalism to systems moving in 
the reference frame considered, but also in order to rec- 
ognize the relativistic covariant structures of |T]) and ^ , 
we have to go over to a four-dimensional (4D) formalism. 
We collect D and H in the 4D excitation tensor density 
©^(L>, H) = and E and B in the 4D field strength 

tensor F^(E, B) = —F vll , with /i, v, ■ ■ ■ = 0, 1, 2, 3 and 
coordinates (x° = t, 3 ), see Post |3j. Then the 

Maxwell equations read 

d^ u = r , d„F uX + d„F^ + d x F^ = , (3) 

with .j^ (p, j) as4-current. This is the "premetric" form of 
Maxwell's equations. They are covariant under general 
coordinate transformations and do not depend on the 
metric of spacetime, that is, they are valid in this form 
in special and in general relativity alike. 

In order to complete the Maxwell equations to a 
predictive physical system, we have to specify a constitu- 
tive law linking <& tiv to F^ v . In vacuum, here now even- 
tually the metric g^ of spacetime enters with signature 
( — h + +), we have 



(4) 



with g := dot g pa 7^ 0. From the covariant components 
of the metric g^ v , its contravariant components g Xv can 
be determined via g^\ g Xl/ — 8" . 



For magnetoelectric media that are local and linear, 
we assume, following Tamm [TtJ and Post Q, the consti- 
tutive law 



^ ~ 2 ' 



(5) 



where x Xl " JK is a constitutive tensor density of rank 4 and 
weight +1, with the dimension [x] — 1/ resistance. Since 
both <5 Al/ and F aK are antisymmetric in their indices, we 



have x = _^AiyK<T _ _^,;yA<TK ^ n an tisymmetric pair 
of indices corresponds, in 4D, to six independent compo- 
nents. Thus, the constitutive tensor can be considered as 
a 6 x 6 matrix with 36 independent components, see (fT|) 
and ©. 

A 6 x 6 matrix can be decomposed in its tracefree 
symmetric part (20 independent components), its anti- 
symmetric part (15 components), and its trace (1 com- 
ponent). On the level of x Aly<TK , this decomposition is 
reflected in 

36 = 20 © 15 © 1. 

The third part, the axion part, is totally antisymmet- 
ric and as such proportional to the Levi-Civita symbol, 
(3) x a^ ctk . = x [Xuctk] = ae XvaK . Here, the totally an- 
tisymmetric Levi-Civita symbol is e\ vaK — ±1,0; we 
denote pseudotensors with a tilde. The second part, 
the skewon part, is defined according to i 2 )-^ 1 '^ 9 := 
2"(x MI/Ap — X ApA "0- If the constitutive equation can be de- 
rived from a Lagrangian, which is the case as long as only 
reversible processes are considered, then ( 2 )^ Al/<T ' i = 0. 
We will assume this condition henceforth. Below, for 
Ci"203, it will be verified experimentally. The principal 



part ('V'" 7 ' 



has the symmetries ( 1 )^ ;Wk 



and W x I^°k] = Q. The tensor x XvaK has now 20 + 1 
independent components and the constitutive law reads 



C5 



1 



( V w,i + ae XvaK ) F aK . 



(7) 



We can express the axion piece a directly in the consti- 
tutive tensor ©. With ( 2 ) x Al «™ = 0, we find 

~ _ 1 ~ , ,\vok _ 1 /,,0123 1 ^0231 , ^0312\ io\ 

4! _ 3^ ' ' 

This 4D pseudoscalar (or axion piece) of ^- Wft will be 
determined for Cr2C>3. 

We split (O into space and time [HI • Then we recover 
equations of the form of ((T|) and ^j, but with an exact 
relativistic meaning (a, 6, • • • = 1, 2, 3): 

D a = e ab E b + 7% B b + aB a , (9) 
H a = [ i a lB h - ^ a E h -aE a . (10) 

We have the 6 permittivities e ab = e ba , the 6 permeabil- 
ities fi a b = f'ba, and the 8+1 magnetoelectric pieces j a t 
(its trace vanishes, 7 C C = 0) and a, respectively. Equiva- 
lent constitutive relations were formulated by Serdyukov 
et al. [IH, p. 86, and studied in quite some detail. It is 
remarkable, as can be recognized from ([7]) and ©, that 
a is a 4D pseudoscalar. 
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FIG. 1: The MEe effect (linear magnetoelectric effect with 
electric field-induced magnetization) of Cr 2 0a: Temperature 
dependence of the magnetoelectric components a\\ and a± 
according to Astrov (14 ] . 



IV. CHROMIUM OXIDE Cr 2 3 AND ITS 
CONSTITUTIVE LAW 

On the basis of neutron scattering data and suscepti- 
bility measurements of the antiferromagnetic chromium 
sesquioxide Cr203, Dzyaloshinskii [l2j was able to estab- 
lish the magnetic symmetry class 3'to' of a Cr203 crystal. 
On this basis, he formulated the following constitutive 
law for Cr 2 03: 



D z 
B z 



£\\£oE z + — H z , 
c 

an 

Mil tJ-oH z H E z . 

c 



(11) 
(12) 
(13) 
(14) 



The z-axis is parallel to the trigonal (and the optical) axis 
of the crystal. The permittivities parallel and perpendic- 
ular to the z-axis are denoted by £||, s±, analogously the 
permeabilities by fi\\,n±, and the magnetoelectric moduli 
by a||,aj_. Note that all these moduli are dimcnsionlcss 
(in all systems of units) . 

The theory pip to (|14p and also the corresponding 
measurements were made in the (D, B) system. However, 
in order to get to the manifestly relativistic covariant 
(D,H) representation ([5]). ([TP]), we have to resolve the 
(D, B) system with respect to D and H and to compare 
the corresponding coefficients. We find magnetoelectric 
matrix 



3 VM± Mil/ U 0-2; 



and the pseudoscalar or axion piece 



(15) 
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FIG. 2: The MEh effect (linear magnetoelectric effect with 
magnetic field-induced polarization) of Cr2 03: Temperature 
dependence of the magnetoelectric moduli an and a± accord- 
ing to Wiegelmann et al. [2(j, Fig.2. Their measured values 
are denoted by large empty triangles, pointing up or down; 
superimposed are the interpolated and digitalized point, small 
black triangles, pointing up or down. The curve was deter- 
mined by using B-splines. One should also compare Rivera 
[2lT | and Wiegelmann [22] . The relative sign between a|| and 
a± was taken from Astrov [lj. The magnetic field fj,oH was 
below 6 tesla. 



3 V Mi- M|| 



(16) 



Eqs. lfToj) and (floj) can be collected in the "relativistic" 
a-matrix 



ot b ■= 7 b + aOj 




(17) 



MEASUREMENTS ON Cr 2 3 



Astrov [14|, see Fig.l, measured an and a± in an 



electric field E according to ([13)) and (| 14[) and Rado & 
Folen in a magnetic excitation H according to (|lip and 
(fT2"[) . Within the measurement limits, they found the 
same values. This verifies that ( 2 )^- Wk = and con- 
firms Dzyaloshinskii's theory (below the spin-flop phase). 
Later on, mostly measurements in magnetic fields below 
UnH = 6T were made, see Rivera [2l[ and Wiegelmann 
et al. 0. 

Measurements of Wiegelmann et al. [2(| are plotted in 
Fig.2. The maximum of ct\\ was found at about 275 K: 



a N (at 275 K) w 4.13— x c 
m 



1.238 x 10" 



(18) 



Now we can compute from the values of Fig.2 the pseu- 
doscalar (|16p . However, we need additionally the perme- 
abilities fin and They have been measured by Foner 
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We find at 4.2 K, [i± « 1.00147 and the maximum 
value near the Neel temperature » 1.00162. Since /in 
deviates even less from 1, we have Then (|16p can 

be easily evaluated: a » i (2ai + an) lo- Our results 
are plotted in Fig. 3. 

As we can see, for temperatures of up to about 163 
K, the pseudoscalar is negative, for higher temperatures 
positive until it vanishes at the Neel temperature of about 
308 K . For the maximum, we find 

a max (at 285 if) w 3.10 x 10" 4 Y » 0.822 . (19) 

We conclude that 5 is fairly small but, for T 7^ 163 if, 
definitely nonvanishing. 
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FIG. 3: Our new result: The pseudoscalar or axion piece 5 
of the constitutive tensor ^ Xl " TK of Cr2 03 in units of Yn as a 
function of the temperature T in kelvin; here Yn is the vacuum 
admittance, which is 1 in Gaussian units and about 1/(377 
ohm) in SI. In the figure the tilde of a is missing. 

Does this result imply consequences also for the ex- 
perimentalist? We think so for the following reason: As 
we saw above, as well as are approximately one. 
Therefore the magnetoelectric 7 matrix (|15p becomes 

y A 0\ 

7°6«-#(o.L-«||) 1 . (20) 
6 \0 0-2 J 

The question is now: Can we find a substance in which 

a_L=Q!||, (21) 

that is, in which the matrix 7% vanishes for all tem- 
peratures? This challenge for experimentalists would be 
interesting in the sense that then one would have a sub- 
stance in which the only magnetoelectric piece would be 
the pseudoscalar (or axion) piece a. In other words, this 
substance would display an isotropic magnetoelectric ef- 
fect. 



For a theoretician it could be of value if he/she looks 
for a microscopic Hamiltonian. The pseudoscalar a of 
the magnetoelectric effect should have a different physi- 
cal origin as compared to the 7 part. Thus, it could be 
helpful for developing microscopic models for the magne- 
toelectric effect. 



VI. PROPERTIES OF THE PSEUDOSCALAR 
OR AXION PIECE 

Unlike the 3D vectors of the electric an the magnetic 
fields and the 3D tensors of permittivity e ab , of perme- 
ability /i ab , and of the magnetoelectric moduli j a b, which 
all depend on the choice of the reference frame and the 
local coordinates, the value of a is always the same. It 
is invariant under any orientation-preserving transforma- 
tion of frames and coordinates — and changes sign when 
the orientation is changed. 

If we consider the pseudoscalar a alone, then we can 
take its constitutive law from © and (fTTJjl . 

D a = +aB a , H a = -aE a (22) 

or, in 4D, 

& xu _ s^Wp^ j 2 . (23) 
A space reflection 

D a - -D\ H a -> H a , E a -> -E a , B a ^B\ (24) 
as well and a time inversion 

D a - D\ H a - -H a , E a - E a , B a - -B a , (25) 

see Janner [24| and Marmo et al. [25( , will turn (|2"2"f into 
its negative, 

D a = -a B a , H a = +a E a . (26) 

This is an expression of the pseudoscaXai nature of a. 
Therefore a is P odd and T odd. 

Moreover, the energy-momentum tensor for the elec- 
tromagnetic field, see [3, HI] , 

1x v = \ ® aT F aT 5l - <5™F Xa , (27) 

if (|2^| is substituted, vanishes: 

1\ v (of axion piece a) = . (28) 

Thus, the electromagnetic energy density To , the 
energy-flux density (Poynting flux) 1 b etc. of the axion 
piece vanish. 
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VII. ANALOGUES OF THE 4D 
PSEUDOSCALAR a IN NETWORK THEORY, IN 
ELECTRICAL ENGINEERING, AND IN 
PARTICLE PHYSICS 

The structure of the constitutive law (f2"2")l or (f2"3")) is not 
unprecedented. In electrical engineering, in the theory 
linear networks, more specifically in the theory of two 
ports (or four poles), Tellegen [1, [26| came up with the 
new structure of a gyrator, which is defined via 

t>i = -si 2 , v 2 = six, (29) 

where v are voltages and i currents of the ports 1 and 2, 
respectively. Let us quote from Tellegen |26|. p. 189: "The 
ideal gyrator has the property of 'gyrating' a current into 
a voltage, and vice versa. The coefficient s, which has the 
dimension of a resistance, we call the gyration resistance; 
1/s we call the gyration conductance." The gyrator is a 
nonreciprocal network element. 

If we turn to the electromagnetic field, then because 
of dimensional reasons the quantities related to the cur- 
rents ii, ii are the excitations D a , H a and the quantities 
related to the voltages Vi, v 2 the field strengths E a , B a . 
Then we find straightforwardly the analogous relations 

E a = -sH a , B a = sD a . (30) 

If we rename the resistance s according to s = 1/5, then 
([30]) and ([22]) coincide. Without the least doubt, the 
gyrator is in the theory of electrical networks what the 
axion piece is in magnetoelectricity. The axion piece 'ro- 
tates' the excitations, modulo an admittance, into the 
field strengths, as the gyrator the currents into voltages. 

These analogies or rather isomorphisms carry even fur- 
ther. In 2005, Lindell & Sihvola Q, see also [23], intro- 
duced the new concept of a perfect electromagnetic con- 
ductor (PEMC). It obeys the constitutive law (|2"2")l or 
(|23| . The PEMC is a generalization of the perfect elec- 
tric and the perfect magnetic conductor. In this sense, 
it is the 'ideal' electromagnetic conductor that can be 
hopefully built by means of a suitable metamaterial, see 
Sihvola [28|. 

Continuing with our search for isomorphisms, we turn 
to axion electrodynamics, see Ni [i"o| and Wilczek [Tlj 
and, for more recent work, Itin [2!|, [30|. If for vacuum 
electrodynamics we add to the usual Maxwell-Lorentz 



expression specified in ([4]) an axion piece patterned after 
the last term in (J7|), then we have the constitutive law 
for axion electrodynamics, 

= Y V=gF Xu + i 5? W F CTK . (31) 

In Cr2C>3 we have a w 10 — 4 Yq. It is everybody's guess 
what it could be for the physical vacuum. In elementary 
particle theory one adds in the corresponding Lagrangian 
also kinetic terms of the axion a la ~ g^d^ad^a and 
possibly a massive term ~ m|S 2 . However, this hypo- 
thetical P odd and T odd particle has not been found 
so far, in spite of considerable experimental efforts, see 
Davis et al. [3lj . 

The axion shares its P odd and T odd properties with 
the a piece of Cr2C>3, with the gyrator, and with the 
PEMC. One may speculate whether an axion detector 
made of Cr 2 03 crystals could enhance the probability of 
finding axions. 



VIII. CONCLUSIONS 

Our results can be summed up as follows: 

• The magnetoelectric pseudoscalar of Cr203 is tem- 
perature dependent and of the order of 10 _4 Y"o- Thus, 
the Post constraint is invalid in general. 

• Our result and the formalism is of general interest to 
areas such as multiferroics and similar materials where 
there is a need for a generalized and concise description 
of the phenomena. 

• We suggest to the experimentalist to search for a 
substance with an isotropic magnetoelectric effect. 

• Since the gyrator and the PEMC are established 
notions in network theory and in electrical engineering, 
these two notions, together with the isomorphic struc- 
ture of the pseudoscalar piece in Cr203, support the ex- 
istence of a fundamental axion field coupled to conven- 
tional vacuum electrodynamics. Thus, axion electrody- 
namics gains plausibility by our results. 
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